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O ■ Abstract 

t> : 

^ Radiative decays of the 4> meson have been studied using a data sample 

of about 19 million cf> decays collected by the CMD-2 detector at VEPP-2M 
collider in Novosibirsk. From selected e + e" - -> 7T 7r°7 and e + e — > ?]7r 7 events 
the following model independent results have been obtained: 

Br((/) 7r°vr 7) = (0.92 ± 0.08 ± 0.06) x 10~ 4 for M^o > 700 MeV, 
^ Br((f> -► r/7r°7) = (0.90 ± 0.24 ± 0.10) x 10~ 4 . 

It is shown that the intermediate mechanism Jq (980)7 dominates in the <p — > 
7r°7r°7 decay and the corresponding branching ratio is 

Br((f> -» / (980)7) = (2.90 ± 0.21 ± 1.54) x 10" 4 . 
The systematic error is dominated by the possible model uncertainty. 

Using the same data sample the upper limit has been obtained for the P- and 
CP-violating decay of rj at 90% CL: 
Br(rj -> vr 7r°) < 4.3 x 10~ 4 . 



X 



Introduction 

In the preceding paper [ [ij we reported on the first observation of the — > tt + tt~j 
decay with the CMD-2 detector. The present work provides complementary informa- 
tion on the electric dipole radiative transitions of the meson obtained from the 7r°7r°7 
and 777r°7 final states. 

The purely neutral decay <fi — > 7r°7r°7 has no bremsstrahlung background as the 
cf) — > 7r + 7r~7 decay and is the most efficient to study the two pion mass spectrum in the 
4> — > /o (980)7 decay. With the CMD-2 detector this mode as well as another decay 
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with five photons in the final state — > i]7c 0, y have been studied and the first results [ 
|2| based on about 25% of the data confirmed those reported earlier by SND [ [3], 

The CMD-2 detector described in detail elsewhere [ [| has been taking data since 
1992. In addition to the barrel Csl calorimeter, the endcap calorimeter made of BGO 
crystals was installed in 1995 making the detector almost hermetic to the photons. 
The energy resolution for photons in the Csl calorimeter is about 8% independent of 
the energy and <Je/E = 4.6%/ \J E(GeV) for the BGO calorimeter. 

In total, the 14.2 pb~ x of data have been collected at 14 energy points around the 
mass. For the analysis of the — > 7r°7r°7 and — > rjir ^ decays presented here 12.8 
pb~ x were used corresponding to 18.8 x 10 6 decays. 

4> — ► 7r°7r°7 channel 
Selection of 7r°7r°7 Events 

Candidates for this decay were selected from a sample of purely neutral events with 
the following criteria: 

1. There are five or more photons in the Csl and BGO calorimeters with the total 
energy deposition E tot > 1.75 • Eb ea m- The Monte Carlo simulation of the 7r°7r°7 
events well reproduces the resolution over E tot and shows that 97% of signal 
events survive this cut. A requirement of minimum three photons in the Csl 
calorimeter gives a high trigger efficiency [ |6| . 

2. All photons have polar angles in the range 0.6-2.54 radians and have energy 
higher than 20 MeV. 

3. To select 7r°7r°7 events a constrained fit requiring energy- momentum conservation 
was performed finding two best combinations of photon pairs with 7r° masses. 
The x 2 distribution after the constrained fit is shown in Fig. [TJa. Events with 
X 2 < 6 were taken for further analysis. This cut suppresses to the 10~ 3 level the 
background from the — > K^K^ decay mode when Kg decays to 7r°7r° and K® L 
produces the fifth cluster in the calorimeter. 

4. Photons from reconstructed 7r°'s in the BGO calorimeter have energy higher than 
40 MeV. This cut removed incorrectly reconstructed events with photons from 
the beam background. 

5. A selection cut | (P^ —P 2 o)/ (-Pjo +P%o)\ < 0.8, where P^ , P 2 o were pion momenta, 
was applied to remove incorrect combinations in which a free photon of low energy 
was used as a part of a reconstructed 7r°. 
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Figure 1: Study of the <fi — > 7r°7r°7 events: a. \ 2 distribution for a five photon sample 
after the constrained fit. Histograms present the normalized simulated contribution 
from 777 events (lower) and a sum of this distribution with that for ujtt events (upper); 
b. 7r°7T invariant mass for events with 2.4 < % 2 < 6. The histogram is simulation of 
777 events ; c. 7r°7 invariant mass. Histograms are as in "a"; d. 7r°7r° invariant mass 
for events with \ 2 < 2.4. Histograms show a contribution from un° events after the 
"anti-cu" cut (dashed) and a sum of this contribution with 777 background (dotted). 

Background Subtraction 

The main background to events of interest comes from the processes — > 777, i] — > 
7r 7r°7r° with two lost photons and from e + e~ — ► um° with the uj — ► tt 0/ j decay. 

The subtraction of the background from 777 events was performed according to 
the simulation. The invariant mass of the 7r°7r° system for 2.4 < % 2 < 6 is shown 
in Fig. [Tjt) and demonstrates good agreement of the observed background spectrum 
with simulation. The normalization of observed background events to the collected 
integrated luminosity gives A^, = (20.6 ± 1.0) millions of decays in agreement with 
the (18.8 ± 0.9) millions obtained from the analysis of seven photon events from the 
~~ > VI ~^ 7r°7r 7r 7 decay [ [/J. The latter number was used for the normalization of 
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signal events. 

The invariant mass of the 7T 7 system for selected events is shown in Fig. |l]c and 
demonstrates the presence of events from e + e _ — ► ujtt process. The number of ujtt 
events was found to be 506 ± 28 with the 7r°7 mass of (782.0 ± 0.9) MeV close to the 
world average value of the u meson mass [ These events were used to check the 
above cut efficiencies and event selection criteria. By looking for the ujtt signal in 
a sample of six and seven photon events it was confirmed that 92% of signal events 
had exactly five photons. This was used for the efficiency correction. The simulated 
detection efficiency for e + e~ — > ujtt events was found to be (14.4 ± 0.1)% . 
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Figure 2: The cross section vs. energy 
for selected ujtt events. 



Figure 3: The simulated efficiency vs. 
7r°7r° invariant mass. The line is an ap- 
proximation used for calculations. 



The energy dependence of the cross section of the process e + e~ — > ujtt — > 7r°7r°7 
is presented in Fig. |2|. It was fit taking into account the interference of the non- 
resonant process with the <fi decay. The values of the parameters of the — > ujtt decay 
were taken from [ [|. The non-resonant cross section obtained at 2Eb eam = is 
o"o = (0.58 ± 0.02 ± 0.04) nb consistent with that measured by SND [ [|. The second 
error represents the systematic uncertainty caused by background subtraction. 

The "anti-o;" cut Mj n „(7r°7) < 750 MeV reduces the admixture of ujtt events to 
the level of about 5% (see the dashed histogram in Fig. |l]d). However, it does not 
completely remove ujtt events mostly because of some incorrectly reconstructed events 
in which a free photon is used as a part of a reconstructed tt°. These incorrect com- 
binations were studied using simulated ujtt events and experimental events from the 
"off-0" region where the process e + e~ — ► ujtt dominates. The ratio of the number of 
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incorrectly reconstructed events to the total number of ujtt events was found to be 
0.060±0.003 for simulation and 0.09±0.02 for experimental data. The difference was 
used to estimate a systematic error so that the final admixture of ujtc events in the 
7r°7r°7 sample was (5. Oil. 5)%. 

Branching Ratio Calculation 

The invariant mass distribution of the 7r°7r° system for \ 2 < 2.4 is shown in Fig. |l]d with 
the expected backgrounds from ujtt after the "anti-u;" cut and from 777 events. The 
distribution demonstrates the increase in the number of events with a high invariant 
mass (a free photon of low energy). In total, 268 ± 27 of 7r°7r°7 events have been found 
after all cuts and background subtraction. 

The detection efficiency as a function of the 7r°7r° invariant mass was obtained 
using simulation of the process e + e~ — > X(M) r ) where X(M) was a particle with a 
small width and variable mass M decaying into 7r°7r°. 

The angular distribution of the free photon was taken as dN/dO^ ~ (1 + cos 2 (# 7 )). 
The detection efficiency obtained is presented in Fig. ^| as a function of the 7r°7r° 
invariant mass. 

The 7r°7r° mass spectrum was obtained from the experimental distribution (Fig. [l|d) 
after background subtraction and taking into account the detection efficiency for each 
histogram bin. The spectrum was normalized to the number of decays obtained 
from the — > 777 — > 7r 7r°7r°7 analysis [ [7| . The resulting differential cross section vs. 
invariant mass is presented in Fig. |]a showing a resonance increase at high masses. 
Figure [|b presents the angular distribution for free photons for signal events with 
M^o^o > 800 MeV. The line shows the distribution expected for a scalar intermediate 
resonance: dN/d6 1 « (1 + cos 2 (6 l 7 )). Open points show the angular distribution of 
the subtracted background. The branching ratio calculated from the integral over the 
differential spectrum was found to be: 

Br((j) -> 7r°7r°7) = (1.08 ± 0.17 ± 0.09) x 10~ 4 
for the whole invariant mass range. The systematic error comes from the uncertainty 
of the background subtraction (about 5%) and from the uncertainty of the number of 
mesons (about 5%). The main contribution to the statistical error comes from the 
region M^o < 550 MeV because of the big uncertainty of the background subtraction. 

This result can also be presented as: 

Br((j) -> 7r°7r°7) = (1.06 ± 0.09 ± 0.06) x 10" 4 for M^o > 550 MeV; 
Br((f) -> 7r°7r°7) = (0.92 ± 0.08 ± 0.06) x 10~ 4 for M^o > 700 MeV; 
Br((p -> 7r°7r°7) = (0.57 ± 0.06 ± 0.04) x 10~ 4 for M^o^o > 900 MeV. 
The results above are consistent with those presented by SND [§,§]]. 

It should be mentioned that only results for M^o > 700 MeV can be considered 
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Figure 4: a. Differential cross section vs invariant mass. The dotted line is the four 
quark model fit with the branching ratio of 3.11 x 10~ 4 . The dashed line is the narrow 
pole fit; b. Angular distribution for signal events and background (open points). The 
line is dN/d9 y « (1 + cos 2 (6> 7 )) 



as model independent. For lower masses the interference with the nonresonant (um° — > 
7r°7r° 7) and resonant (0 — > /wr — > 7r°7r 7) backgrounds can change the branching ratio. 
To extract these contributions a much better measurement at lower masses is needed. 

About half of the signal events have the free photon energy below 100 MeV pointing 
to the presence of the /o(980) intermediate state in the 7r°7r° system. 



Data Interpretation 

For data interpretation, similarly to our analysis of the <p 7r + 7r _ 7 decay [ the 
model suggested in [ [II], [□]] was used. The model considers the <fi — > fo (980)7 decay 
under different assumptions about the f meson structure (a two-quark or four-quark 
state or KK molecule depending on the values of the coupling constants g^/47r, 
g 2 T /47r) and calculates the differential cross section over the invariant mass of two 
pions M nx . Also considered are effects of other possible intermediate states like and 
p7r° — > 7r°7r°7 with the total contribution of about 15%. 

The experimental spectrum Fig. f|a was fit using the differential cross section from [ 



10 , 11 1 assuming that the — > decay is completely dominated by the fo (980)7 



mechanism. The following model parameters have been obtained: m/ = (977 ± 3 ± 
6) MeV, g^/47r = (2.44 ± 0.73) GeV 2 , g 2 7r /4 7 r = (0.40 ± 0.06) GeV 2 , Br((j) -> 
/ (980)7) = (3.05 ± 0.25 ± 0.72) ■ 10" 4 . 

The second error includes experimental uncertainties at high 7r°7r° masses and the 15% 
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effect of other processes. The parameters obtained are in good agreement with those 
presented by the SND group [ f| and can be compared to the values obtained from our 
analysis of the <p -> 7r+7r~7 channel [|: m fo = (976 ±5 ±6) MeV, Br{(f> -> / (980)7) = 
(1.93 ±0.46 ±0.50) • 10" 4 . 

The results for the m/ and Br(<p — > / (980)7) obtained from the 7r°7r°7 and 7r + 7r~7 
final states are consistent within the errors. Therefore, we can perform their combined 



analysis. A simultaneous fit of the photon spectra using the formulae from [ [LI], [12], P~3 
allowed to determine the model parameters with a better statistical accuracy: 
m fo = (975±4±6) MeV, g^/4vr = (1.48±0.32) GeV 2 , g^/47r = (0.41±0.06) GeV 2 , 
Br((p -> / (980)7) = (2.90 ± 0.21 ± 0.65) ■ 10~ 4 , and # = (1.47 ± 0.19) radians. The 
X 2 /d.f.=1.5 has been obtained. The results of the fit are shown in Fig. ^a by the 
dotted line. According to the model, the obtained values of the coupling constants 
(or the branching ratio) could only be explained if /o is a four quark state. The 
ratio R=g 2 KK /gl n had weak dependence on the f structure and was found to be 
R=3.61±0.62. 

In the above interpretation the whole visible signal was due to the <fi — > /o7 decay 
whereas the influence of other possible mechanisms (pn, 0-7) was estimated to be about 
15% and included into the systematic error. 

To understand the sensitivity of the branching ratio of the — > fo (980)7 decay to 
the particular model, we used another approach and performed the narrow pole fit. 
The 7r°7r° mass spectrum was fit with the function: 



da 911/ f I Mgg \ ft 4m l I m / (980) r / (980) I ^e"^m /(1200) r /(1200) i 2 

dM ^ MV1 ™ \ l i^i-y 1 \ A /(980) + A /(1200) I ' 

where A/ = M 2 ^— mj + iM 7rn T(M 7rn ). The parameters of the /o(1200) (or a) resonance 
could not be extracted from our data and were fixed at mj( 12 oo) = 1200 MeV and 
r./(i2oo) — 600 MeV. These parameters have very small influence on the result and 
instead of the /o(1200) Breit-Wigner the constant background amplitude Ae 1 ^ can be 
used. The relative phase ip was found to be close to zero and fixed at that value. 

The following parameters have been obtained: 
mfo = (987 ± 7) MeV/c 2 , T /o = (56 ± 20) MeV, the relative amplitude A=0.22±0.09 
(or A=0.08±0.03 in case of constant background), and x 2 /d.f.=0.96. The result of the 
fit is shown in Fig. [|a by the dashed-dotted line. 

The corresponding branching ratio (assuming that only 1/3 is seen in the 7r°7r° 
mode) Br(<p -»■ / (980)7) = (1.5 ± 0.5) ■ 10" 4 has been found. 

The discussion above shows that the data obtained cannot be interpreted without a 
resonance in the two pion mass spectrum at about 980 MeV both in charged and neutral 
modes and the <fr fol branching ratio cannot be lower than about 1.0- 10~ 4 . As many 
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authors agreed [ H , a value that high can hardly be explained in the frame of the two 
quark model of the /o(980) which predicts the branching ratio at the level of 0.5 ■ 1CT 4 . 
To study the influence of other intermediate mechanisms, better measurements of the 
7T7T mass spectra are needed. 



Search for <fi — > rjir 7 Decay 

Using the five photon final state, one can also look for the decay mode — > rjir '-/ 
appearing when 77 decays into two photons. This analysis was based on the same 
event sample as in the study of the 7r°7r°7 decay mode with the selection criteria 
described above. 

A fit finding the best combination of photon pairs with the ir° mass constrained 
and reconstructed pion momentum less than 350 MeV/c was used. An additional 
requirement is that the difference in the energy of the photons from the found 7r° is less 
than 80%. It rejects background from low energy photons and increases the efficiency 
for events with a low energy free photon because of a smaller probability of wrong 
combinations. 

The invariant masses of the remaining most energetic photon pairs are shown in 
Fig. §a. 

The main source of background to observed events is the decay — > 777 — > 37r°7 
with two photons lost. The cut 760 < M w o^ < 805 MeV from the reconstruction 
almost removes the background from ujtt . The ujtt events remaining after the "anti-o;" 
cut are shown by the dotted histogram. 

Over the broad background distribution the excess of 80 ± 22 events is observed 
at (545 ± 4) MeV compatible with the mass of the 77 meson. The solid histogram 
shows the sum of the simulated background from — > 777 — > 37r°7 and remaining ujtt 
events normalized to the number of decays. These simulated events were used for 
background subtraction. The dashed histogram at 550 MeV shows a simulated signal 
from the — > f]TT 0/ ~f decay at the 1 x 10~ 4 level. 

The distributions over cosO^ and invariant mass M m ° are shown in Figs. [5|b,c for 
events with 510 < M 77 < 590 MeV after background subtraction and taking into ac- 
count the detection efficiency obtained by simulation (Fig. |5]d) . The following branching 
ratio: 

Br(4> -> 77vr°7) = (0.90 ± 0.24 ± 0.10) x 10~ 4 

has been obtained. The systematic error comes from the uncertainty of the background 
subtraction and of the number of 0's taken for normalization. 

To check the stability of the result these distributions were also obtained in the 
case of "strong" cuts when the additional requirement for the reconstructed photon 
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Figure 5: The (f) — > r]7r 0/ j study. Open points correspond to the "strong" cuts: a. In- 
variant masses of two photons with the highest energy for data (points) and simulation 
(histograms); b. cos9 7 distribution for events around the rj mass after background sub- 
traction. The line corresponds to the expected distribution dN/dO^ ps (1 + cos 2 ^)); 
c. rjn mass distribution; d. Simulated detection efficiency vs. invariant mass. 



with the highest energy was applied: E imax / 'Eb eam < 0.75. This cut reduced the main 
background from — > 777 — > 37r°7 by a factor of 4, but at the same time considerably 
reduced the detection efficiency at higher r)7r° invariant masses (open points in Fig. |5|d). 
The number of observed signal events dropped to 37 ± 12 and the obtained branching 
ratio was the same within statistical errors. 

The obtained invariant mass distribution shows the growth of the cross section to 
higher masses supporting the hypothesis about the ao(980) intermediate state. The ob- 
tained value of the branching ratio can be explained in the four quark model suggested 
in [0,0. 
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Search for 77 — ► tttt Decays 



The selected 7r°7r°7 events can be used to search for the P- and CP- violating decay 
i] — ► 7r°7r°, where the rj comes from the radiative <fi — *■ 777 decay. From 18.8 million </> 
decays used for the analysis of the 7r°7r°7 channel one could expect about 236,000 events 
decayed via the 777 channel. Such a P- and CP-violating decay should be observed as a 
peak in the invariant mass of two pions at M n7T = m v . Figure |] shows the experimental 
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Figure 6: Search for 77 — > 7r°7r° decay. The histogram is simulation. The line shows a 
possible signal at 90% CL. 

7r°7r° mass distribution from selected 7r°7r°7 events. The line corresponds to a fit with 
a linear function and gaussian distribution representing a possible signal at 90% CL. It 
was found that a possible signal does not exceed 11 events for neutral decay mode. The 
histogram shows a simulated signal from the 77 — > 7r°7r° decay at 90% CL. The detection 
efficiency found by simulation was 0.108. The following result has been obtained: 

Br(r/ -> ttV) < 4.3 x 10" 4 , 
which should be compared to the best previous limit 7 x 10~ 4 for the neutral decay 
mode [ . 



Conclusions 

Using 12.8 pb _1 of data collected around the <j) meson (about 19 million decays) 
events with five photons were selected. The reconstruction of the 7r°7r°7 and ry7r°7 final 
states gives the following model independent results: 

Br((j) -> 7r°7r°7) = (0.92 ± 0.08 ± 0.06) x 10" 4 for M^o^o > 700 MeV, 
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Br((j) -> r/7r°7) = (0.90 ± 0.24 ± 0.10) x 10" 4 , 
in agreement with those reported by the SND group [ [|, |3| . A resonance at 980 MeV 
is observed in the 7r°7r° mass spectrum. 

For the 7r°7r°7 channel the branching ratio obtained by integration over the whole 
mass spectrum was found to be 

Br((j) -> 7r°7r°7) = (1.08 ± 0.17 ± 0.09) x 10" 4 , 
but this result can be affected by the possible interference with um Q — > 7r°7r°7 or 
p7r° — > 7r°7r°7 intermediate states at M^o^o < 700 MeV. 

In the 7T7T7 channel the — > /o (980)7 mechanism dominates and the most consistent 
description for both charged and neutral modes can be obtained in the four quark model 
with the branching ratio 

Br((p -> / (980)7) = (2.90 ± 0.21 ± 0.65) x lO" 4 . 
However, this result is model dependent and it is shown that in another approach the 
branching ratio can decrease to about 1.5 x 10~ 4 , still higher than the prediction of the 
two quark model. The difference between two approaches can be used as an estimate 
of the additional systematic uncertainty because of the model so that 

Br((p -> / (980)7) = (2.90 ± 0.21 ± 1.54) x 10" 4 . 

From the discussed models the following / (980) parameters can be obtained: 

m fo = (978 ± 4 ± 8) MeV, T fo = (56 ± 20 ± 10) MeV, 
where the second error represents a model dependent uncertainty. 

In the 7]7T 0/ ~f final state events with the high rj-ir invariant mass dominate. This 
supports the hypothesis about the ao(980) intermediate state. The obtained branching 
ratio is also higher than the prediction of the two quark model. 

For the P- and CP-violating decay of rj —>■ 7r°7r° the following upper limit at 90% 
CL has been obtained: 

Br(r) -> 7r°7r°) < 4.3 x 10~ 4 . 
This result is the most stringent upper limit today. 

The values of the branching ratios Br((f> — > / (980)7) and Br((f> — > a (980)7) 
obtained in this work are rather high and according to theoretical predictions favour 
the four-quark structure of these states [ 10, 14 . However, predictions for the branching 
ratios significantly differ from one work to another and can hardly be considered as very 
accurate. Moreover, recently there were successful attempts to explain the observed 
distributions in radiative decays <fi — > tttt^ as well as their branching ratios without 
assumptions about the exotic structure of the / (980) meson [ [Ll| ITT . 



In general, the situation with the identification of the members of the lowest scalar 
nonet and the related problem of the /o(980) and a (980) structure remains to be 
unclear and controversial. While Ref. [ 113 presents many arguments in favour of the 



four-quark structure of these states, other possibilities are widely discussed. Ref. [ [19 
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advocates the idea that the /o(980) meson doesn't belong to the qq family and originates 
from the KK molecule embedded in the continuum. However, in [ 20] the idea of the 
loosely bound KK molecule is strongly disfavoured and the possibility is discussed that 
both / (980) and a (980) are conventional qq states with properties strongly distorted 
by coupling to the nearby KK threshold. The small value of the two-photon width of 
these states which traditionally used to be one of the strong arguments in favour of their 
exotic structure seems to be well accounted for in different theoretical models assuming 
the conventional qq structure of / (980) and a (980) mesons [ [16], ^l], |22| . Note also 



the results of OPAL [ |23| and DELPHI [ ^4| who studied the inclusive production of 
/o(980) and ao(980) mesons in hadronic Z° decays and found that their properties are 
consistent with the normal qq mesons. 

One can conclude that the problem of the / (980) and a (980) structure is far 
from being solved and requires both new theoretical approaches and reliable, model 
independent experimental results. 
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